load('Data')

Q=diag([1000 4000 0 0]);
R=10;
[K,S,e]=lqr(A,B,Q,R);